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^_> ' Abstract 

a • 

In the paper we study the subject of positivity of systems with se- 
quential fractional difference. We give formulas for the unique solutions 
to systems in linear and semi-linear cases. The positivity of systems is 
considered. 



> 

00 ' 1 Introduction 

^- 

^l ' The first definition of the fractional derivative was introduced by Liouville and 

Riemann at the end of the 19-th century. Later on, in the late 1960s, this 
idea was used by engineers for modeling various processes. Thus the fractional 
calculus started to be exploited since that time. This calculus is a field of 
mathematics that grows out of the traditional definitions of calculus integral and 
derivative operators and deals with fractional derivatives and integrals of any 
order. Fractional difference calculus have been investigated by many authors, 
£ , for example, [1,3,4,6,8,10,12,14-17] and others. The subject of positivity is 

Cy_' well developed for fractional linear systems with continuous time, see [11-13]. 

However, positivity of fractional discrete systems with sequential /i-differences 
is still a field to be examined. In the present paper we open our studies in 
this field. We give formulas for the unique solutions to systems in linear and 
semi-linear cases. Moreover, the positivity of systems is considered. 
The paper is organized as follows. In Section 2 all preliminaries definitions, 
facts and notations arc gathered. Section 3 presents systems with sequential 
fractional differences with results on uniqueness of solutions. Semilincar systems 
we included in Section 4. The last Section concerns positivity of considered 
systems. 
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2 Preliminaries 

Let us denote by To the set of real valued functions denned on D. Let h > 
0, a > and put (hN) a ■= {a, a + h, a + 2h, ...} for h > and ael. Let 

Rl = {i£K":j,>0,1<8< n}. 

Due to notations from the time scale theory the operator a : (hN) a — > (hN) a 
is defined by a(t) := t + h. The next definitions of /i-difference operator was 
originally given in [5], but here we propose simpler notation. 

Definition 1. For a function x £ F(hn) a the forward h-difference operator is 
defined as 

/ . w s x(a(t)) — x(t) 

(A h x)(t) := v v ; ; ii t = a + nh,n£fi , 

h 

while the h-difference sum is given by 

n 

(aA^i) (f) :=ft5^x(o + fcft), 

fc=0 

where t = a + (n + l)h, n e N and (aA^x) (a) = 0. 

Definition 2. [5] For arbitrary f,a£R the h-factorial function is defined by 

* r^ + i-a)' Uj 

where T is the Euler gamma function, j- Z_ := { — 1, —2, —3, . . .}, and we use 
the convention that division at a pole yields zero. 

Notice that if we use the general binomial coefficient (^) := pfb+nrf -h+ii ' 
then (1) can be rewritten as 



#°=fc°T(a + !)(*) 



The next definition with another notations was stated in [5]. Here we use 
more suitable summations. 

Definition 3. For a function x G J~(hK\ a the fractional h-sum of order a > is 
given by 

( a A~ a x) (t) := -— J2(t ~ <K« + fcfc))£* _1) *(« + M) , 

where t = a + (a + n)h, n £ N . Moreover we define ( a A°x) (t) := x(t). 
Remark 4. Note that a A^ a : J-(hN) a ~ > -^"(ft,N) a+Qh - 



Accordingly to the definition of /i-factorial function the formula given in 
Definition 3 can be rewritten as: 

(A-.)M=>'E rl xr t ''i) 'i-^) 

fc=0 v 7 

for t = a + (a + n)h, n G No. 

Remark 5. In [10] one can find the following form of the fractional /i-sum of 
order a > 0: 

that can be useful in implementation. 

The following definition one can find in [2] for h = 1. 

Definition 6. Let a G (0, 1]. The Caputo h-difference operator a A^ + x of order 
a for a function x 6 J~<hN) a is defined by 



( a A^x) (t) := ( Q A- (1 - tt) (A h x)) (t), t G (ftN)a+(i-a)fc • 

Remark 7. Note that: „A^ : .F(/,N) a ->■ - 7r (ftN) a+(1 _ Q) , l , where a e (0, 1]. 

We need the power rule formulas in the sequel. Firstly, we easily notice that 
for p^O the well defined /i-factorial functions have the following property: 

A h (t - a)^ = p(t - a)t l K 

More properties of h- factorial functions can be found in [9] . In our consideration 
the crucial role plays the power rule formula presented in [7] , i.e. 

where ip(r) = (r - a + tih)^, r G (hN) a , t G {hN) a+ah . Note that using the 
general binomial coefficient one can write (2) as 

( Q A^)(*) = r(/, + i)( n + " + M )^- 

Then if ij) = 1, then we have for fj, = 0, a = (1 — a)h and t = nh + a + ah 



i-\-a 



r(n + a4 1) _ hQ = ( n + a \ h o 



r(a + l)r(n+l) V n 

Let us define special functions, that we use in the next section to write the 
formula for solutions. 



Definition 8. For a, (3 > we define 



**>*H 0, "" fc j farn<* (3) 



and 



^(n/i := i \ n ' r( M )r(n+i) > ° (4 

10, for ?i < 

where n, fc, s G No and /i = fca + s/3. 
Remark 9. It is worthy to notice that: 

(a) ip 0fi (nh) = 1; 

(b) pi,o(n/0 = C^; 1 )/,- = ( A,; Q 1) ((n - l)/» + a/i); 

(c) ip k ,s{nh) = r ( fcoH^g +i|r (" _ fc + 1 ) anci as tlie division by pole gives zero, the 
formula works also for n < k, n G N; 

(d) ip k ,s( nh ) = r(fca+ s /3+i) ((" " k ) h + kah + #)i fa+S ■ 

We also need the property presented in the following proposition. 
Proposition 10. Let a, (3 G (0, 1], h > and a = (a - l)h, b = ((3 - l)h. Then 

(oA,;>t, s ) (nh + a) = ip k+liS (nh) (5) 

and 

U^h <Pk,8 ) (nh + b) = ip k . s+1 (nh) . (6) 

Proof. We show only equality (5), as (6) is a symmetric one. 

Let n :— ka + sf3. For r <G (hN)kh we define the following n-factorial function 

V'(r) :=(r + /x/i)^ ) . Since 

Pk.-M) = r(fca+ 1 s/?+1) ((« - fc)h + kah + s /3h)i ka+s ^ 
-ip(nh — kh) 



for n > k and (pk tS (mK) = for m < k, by (2) we get 

(oA£<Vm) (*) = K _a ^.) (*) - f^Vl) (oA,rv,) (f) 



1 F ^ + 1 ) (t + uh)^ +a) 



1 



r(/j + a 



T) (t + « 



(m+") 



where t = a — kh + nh. Hence 

t + uh = nh -(k + l)h + (k + l)ah + sj3h 
and 

(oA~> M ) (nh + a) = r0i+ X Q + 1) (nh + a + uh)^ 

T(a + n-{k + l)+ f i + l) _ hll+c 



r(n + a + l)T(n-(k + l) + l) 
T(a + n-k + n) ^ +Q 



T(p + a+l)T(n-k) 

k - 1 + \x 

n — k — 1 



n-k-l + /j, + a\ ^ +a 



fn-(k + l) + (k + 1)q + s /3\ {k+1)a+sfj 

{ n -(k+i) r 

tpk+i, s (nh). 



D 



From the application of the power rule follows the rule for composing two 
fractional /i-sums. The proof for the case h = 1 one can find in [10]. For any 
positive h > we presented the prove in [9]. 

Proposition 11. Let x be a real valued function defined on (hN) a , where a, h g 
R, h > 0. For a, f3 > the following equalities hold: 

( a+ph A- h a ( a A-^)) (t) = ( a A- ( «+^) (t) 

= ( a+ah A-P( a A- a x ))(t), 

where te (hN) a+{a+f})h . 

The next proposition gives a useful identity of transforming Caputo frac- 
tional difference equations into fractional summations for the case when an 
order is from the interval (0, 1]. 

Proposition 12. [9] Let a € (0, 1], h > 0, a = (a — l)h and x be a real valued 
function defined on (hN) a . The following formula holds 

( A^ a (oAg f .x)) (nh + a)= x(nh + a) - x(a), n G Ni . 

3 Systems with sequential fractional differences 

Let a,/3 6 (0, 1] and x : (hN) a -> R". Moreover, let us take a = (a — l)h and 
b = (/? — l)/i. Then we define 

y(nh + b):=( a A^x)(nh). 



Note that y : (hN) b — > 1™. Then we apply the next difference operator of order 
/3 on the new function y and consider here an initial value problem stated by 
the system: 

( a Al,x) (nh) = y(nh + b), (7a) 

( 6 A^y)(n/i) = f(nh,x(nh + a)) (7b) 

with initial values: 

( a A^x)(0) = x 0> (8a) 

x(a) = x a , (8b) 

where x a ,Xo are constant vectors from W 1 . We use Proposition 12 twice and, 
for n > 1, we get: 

(oA- 3 (bAfo)) (nh + b) = y(nh + b) - y(b) 

= („A^,i) (nh) - xq 

and 

(oA" Q ( A^,x)) (n/i + a) = a;(nft + a) - x a . 

Hence 

(aA^x) (nh) =x Q + (oA-^/J (nfc + 6) , (9) 

where f(nh) := f(nh, x(nh + a)). Ncxtly 

z(n/i + a) = x a + x ( A^ a l) (nh + a) + (oA^g) (nh + a) , (10) 

where g(nh) = (oA^/l (nh + b). 

Firstly we prove the formula for the unique solution in linear case, i.e. when 
in equation (7): f(nh,x(nh + a)) = Ax(nh + a), where A is a constant square 
matrix of degree n. 

Theorem 13. The solution to the system 

( a A^x)(nh) = y(nh + b), (11a) 

(bA^yj (nh) = Ax(nh + a) (lib) 

with initial conditions (8), i.e. ( a A^ + xj (0) = xq and x(a) = x a , xo, x a G K™ , 
is given by the following: 

-j-oo 

x(nh + a) = ^2 Ak (Vk,kX a + <Pk+i,kXo) (nh) , (12) 

fc=0 

for n e No . 



Proof. For n = let us notice, that only 930,0(0) = 1, for any k > the next 
terms are zero, so in fact we have: x(0h + a) = x a . 
For n > let us define the following sequence 

x m +i(nh + a) = x a <pofi{nh) + x ip lfi (nh) 

+ ( A- a g m ) (nh + a), m E N , 

where g m (nh) = I oA^ p / m ) (nh + b) and f m (nh) = Ax m (nh + a) with xo{nh + 

a) = x a - 

We calculate the first step. As fo(nh) = Axo(nh + a) = Ax a , then go(nh) = 
Ax a (oA _,3 l) (nh + b) = Ax a tpo,i(nh). Going further, 

xi(nh + a) = x a ipofi{nh) + x (pi t0 (nh) + ( A^ a g ) (nh + a) . 
What could be written as 

xi(nh + a) = (xa<Po,o + x ipi,o + Ax a ipi t i) (nh) . 
and, using Proposition 10, we get 

x 2 (nh + a) = (x a (p Qfi + x <pi, Q + Ax a <Pi,i + Ax Q ip 2 ,i + A 2 x a ip 2 ,2) (nh) . 
and the next element of the considered sequence has the following form: 

x 3 (nh + a) = (x a (po,o + x (pi fi + Ax a <pi,i + Ax tp 2 ,i + A 2 x a if 2 , 2 
+A 2 x (p 3 , 2 + A 3 x a <p 3:3 ) (nh) . 

Taking m tending to +cxd we get formula (12) as the solution of (11) with initial 
conditions (8). □ 

3.1 Semilinear sequential systems 

Firstly we state some technical lemma and notations. 

Lemma 14. Let u : (hN) -4 K and a > 0. Let (oA^^) (nh + kah) = 
71 (nh + kah) and 71 (nh) := 71 (nh + kah) for k 6 N. Then for k £ N we get 

( A^7i) («) = (<A; (fc+1)Q 7 ) (t + kah) , (13) 

where t = nh + ah. 

Proof. First let us consider the case k = 1. Then from Proposition 11 we can 
write 

( ah A- a (0V7)) (t) = (oA- 2 " 7 ) (t), 

where t = nh + 2ah, n £ Nq. 



Let jiinh + ah) = (oA /i Q 7j (nh + ah) and ji(nh) := 71(71/1 + ah). Then 
( Q A- a ^ 1 )(nh + ah) = 



T(a) 
h 



J2 {nh + ah~ a{rh))[ a - 1] ^(rh) 



r-0 
n-\-a 



J2 (nh + 2ah - o-(s/i))£* _1) 71 (sh) 



T(a) ^ 
= ( Q/l A^ Q 7 i) (nh + 2ah) 
= (oA,; 2Q 7 ) (nh + 2ah) . 

The equation (13) for k > 1 follows inductively. □ 

Note that 

(oA- feQ 7 )(n/i + fca/i) 

h " 

= tyTT £ ( nh + kah - °( rh ))h a ~ 1] iW ■ 

Similarly to the procedure presented in the proof of Lemma 14 we can prove 
that for k, s G N and a > 0, /? > 0: 



(oA^-sO ( n h + kah + s(3h) 



li 



T(ka + s(3) 



J2 ((n + ka + s(3-r-l) h^+'P- 1 ) ^ rh ) 



= h^± ^ n - r B l r k ( a + S ^ ] 1 (rh) 
*— ' 1 (ka + spjl [n — r + 1) 

= £ (" " " + n kC l + r SP ~ l ) hka+S "^h) . (14) 

r=0 ^ ' 

Taking u = ka + s/3 and using formula (4) we can write (14) shortly in the 
following way 

n 

(oA^-f) (nh + uh) = ^ <fk,s(nh - rh)j(rh) . (15) 

r=0 

Moreover, we can also write direct formula for values (oA^ a g) (nh + a) given 
in (10) for nonlinear problem. In fact using Definition 3 of fractional summation, 
formula (4) of functions tfk,s and Proposition 11 we write (15) as follows: 

x(nh + a) = x a + x ( A^ a l) (nh + a) + V^ (fii(nh — <r(rh))f(rh, x(rh + a)) . 

r=0 



Using the power rule formula for u = and by Remark 9 we can write the 
recursive formula for the solution to nonlinear problem given by equations (7) 
and conditions (8): 

n 

x(nh + a) = x a + XQ(pi t o(nh) + \. ^Pi,i{nh — o~(rh))f(rh, x(rh + a)) . (16) 

r=0 

The given formula (16) also works for n = as ty?i,i(— h) = 0. Then x(0h + a) = 
x a - We can also check for example the next step: 

x{h + a) = x a + awi.oO) + £i,i(0/i)/(0, x(a)) = x a + x h a + h a+l3 f(Q, x(a)) . 

For special semilinear case when f(nh,x(nh + a)) — Ax(nh + a) + 7(71/1) we 
have f(0,x(a)) = Ax{a) +7(0). Then 

ar(ft + a) = (I + h a+/3 A) x a + h a x + h a+ ^(0) . 

Theorem 15. The solution to the system 

(„A£,s)(nfc) = y(nh + b), (17a) 

( b Af ^j/) (nh) = Ax(nh + a) + -/(nh) (17b) 

wii/i initial conditions (8), i.e. n AjJ t i (0) = xq and x(a) = x a , xo,x a G R" 
is given &y 

x(n/i + a) = ^2 Ak {Vk,kX a + <Pk+i,kXo) (nh) 

k=0 

n / oo \ 

+ E E Ak{ Pk+i,k+i {nh - a(rh)) j(rh) , (18) 

r=0 \k=0 ) 

for n G N . 

Proof. For n = let us notice, that only ipo t o(0) = 1, for any k > the next 
terms are zero, so in fact we have: x(0h + a) = x a . 

For n > based on the proof for linear case we can write the solution formula 
as follows: 



x(nh + a) = ^A k ((p k ,kX a + Vk+i,kXo) (nh) 

fc=0 
-foo 

+ £ A k ( Ai fe+1)(a+ ^ 7 ) (nfc -h + (k + l)(a + /3)/i) 



fc=0 



Then taking into account the formulas (14) and (15) we get the form (18) as 
the solution of (17) with initial conditions (8). □ 



4 Posit ivity 

Based on [12, 13] we consider the following definitions. 

Definition 16. The fractional system (7) is called positive fractional system if 
and only if x[nh + a) e R" for any initial conditions x a , xq € R>. 

Let x , x a e R5 and X n := {X : N ->• W 1 }. The operator T Xo>Xa : X n -)• #„ 
is defined as 

+ OC 



(T^^X) (n) ■= ^2 Ak (Vk,kX a + <Pk+i,kXo) (nh) 

k=0 

n / oo \ 

+ D ( Yl Ak Vk+i,k+x {nh - a{rh)) j(rh) , 

r=0 \fe=0 / 

where X £ X n , X(n) = x(nh + a) and functions tpk,s, (Pk,s are given by (3) and 
(4), respectively. 

The proof of the following proposition is analogous to similar fact in [12]. 

Proposition 17. Let xq, x a S R> and the right hand side of system (7), 
/ fulfils f(nh,x(nh + a)) > 0. ThenT Xo>Xa (R™) C R$. 

Definition 18. The fractional system (7) is called locally positive fractional 
system if and only if for any initial conditions x a , xq S R™ = {x £ R n : Xi > 
0, 1 < i < n} there is r > 1 such that x(nh + a) £ R" for n e [0, r]. 

Proposition 19. J/t/ie matrix I + Ah a+ ^ is positive and 7(71/1) > /or n € No, 
/j/ien /or any xq, x a G R> £«e fractional system (17) is locally positive. 

Proof. Since x(/i + a) = Sfc=o A k (ip k ^ k x a + ip k+1 . k x )(h) = x a + V5i, (/i)x + 
^4(/3i i i(/i)a; a = x a + /i Q a;o + Ah a+ @x a , then with our assumptions, we get local 
positivity. 

□ 
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